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ABSTRACT 

o : 

We present far-ultraviolet spectra of the classical T Tauri star RU Lupi covering the 
912-1710 A spectral range, as observed by the HST/STIS and FUSE satellites. We 
' use these spectra, which are rich in emission and absorption lines, to probe both the 

Q\ ■ accreting and outflowing gas. Absorption in the Lya profile constrains the extinction 

to Ay ~ 0.07 mag, which we confirm with other diagnostics. We estimate a mass 
t-h ■ accretion rate of (5 =L 2) x 10 -8 M & yr _1 using the optical-NUV accretion continuum. 

^f" , The accreting gas is also detected in bright, broad lines of C IV, Si IV, and N V, which 

, all show complex structures across the line profile. Many other emission lines, including 

those of H-2 and Fe II, are pumped by Lya. RU Lupi's spectrum varies significantly 
^) ' in the FUV; our STIS observations occurred when RU Lupi was brighter than several 

other observations in the FUV, possibly due to a high mass accretion rate. 

'"q^! Subject headings: accretion, accretion disks — circumstellar matter — line: iden- 

tification — stars: individual (RU Lupi) — stars: pre-main sequence — ultraviolet: 
$h ' stars 

> ■ 

•r-( , 1. INTRODUCTION 

Classical T Tauri stars (CTTSs) are likely prototypes of the young Sun when it was accreting 
material from its circumstellar disk. In the standard magnetospheric accretion model, the magnetic 
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field truncates the disk at a few stellar radii. Material from the disk accretes onto the star along 
the magnetic field lines and is shocked as it decelerates at the stellar surface. A fast outflow is 
likely associated with the accretion process, and a slow outflow probably emanates from the disk 
(see discussion in, e.g., Hartigan et al. 1995). 

Far-ultraviolet (FUV) spectra reveal these diverse processes. Johns-Krull, Valenti, & Linsky 
(2000) and Calvet et al. (2004) found a correlation between the mass accretion rate and the strength 
of C IV emission from CTTSs. C IV and other high temperature emission lines are much weaker in 
non-accreting naked T Tauri stars (NTTSs) than in CTTSs (Johns-Krull et al. 2000; Ardila et al. 
2002a). Thus, these lines, which trace gas at 10 5 K, are likely produced by the accretion shocks. 

Ardila et al. (2002b) detected wind absorption in the profiles of Mg II lines in HST/GBRS 
spectra of 8 CTTSs, sampling a range of inclinations. They found that the equivalent width of the 
wind absorption is larger for stars with disks that are observed face-on than for stars with disks 
observed edge-on, suggesting some collimation in the outflow. Similar P-Cygni profiles of strong 
lines of O I and C II are typically found in FUV spectra of CTTSs. Fluorescent H2 emission may 
trace both warm molecular gas at or near the disk surface (Herczeg et al. 2002) and molecular gas 
in the surrounding cloud that is shocked as it interacts with stellar outflows (Walter et al. 2003; 
Saucedo et al. 2003). 

In this paper, we present HST/STIS and FUSE spectra, covering the 912-1700 A spectral 
range, of RU Lupi, one of the most heavily veiled 1 CTTSs at optical wavelengths. The spectrum 
is rich in strong absorption features and in emission lines produced by hot and cool gas. In §2, we 
describe the properties of RU Lupi. In §3, we describe our ffi>T/STIS and FUSE observations. In 
§4, we measure the total hydrogen column density in our line of sight to RU Lupi to estimate the 
extinction. In §5, we measure the mass accretion rate and present hot emission in lines that arc 
likely produced by the accreting gas. In §6, we show wind absorption in line profiles. Finally, in §7, 
we analyze variability in the FUV emission from RU Lupi. Our conclusions are presented in §8. 

2. Properties of RU Lupi 

RU Lupi (HD142560) has some of the strongest accretion signatures among CTTSs, including 
strong optical veiling and an Ha equivalent width that has been observed to vary in the range 140- 
216 A (e.g. Appenzeller et al. 1983; Lamzin et al. 1996). Stempels &: Piskunov (2002) estimate an 
optical veiling of 10-20 at 4091 A and 2-5 between 5300-6300 A. This strong veiling, the result of a 
high accretion rate, makes it difficult to estimate stellar properties from the uncertain photospheric 
emission. Stempels & Piskunov (2002) extracted the optical photospheric spectrum of RU Lupi 
from high-resolution VXT/UVES spectra, and measured a radial velocity of — 1.9 ± 1 km s _1 . They 
estimate T e ff = 3950 K and log g = 3.9. Using these numbers, and an approximate mass of M = 0.8 



1 ratio of accretion flux to photospheric flux, see Joy (1949); Herbig (1962) 
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Mq (Lamzin et al. 1996), they calculate R = 1.66 R Q and L = 0.49 L Q , assuming a distance of 200 
pc. The distance is probably closer to 140 pc (de Zeeuw et al. 1999; Bertout et al. 1999). Based on 
our extinction estimate Ay = 0.07 (see §4), d = 140 pc, and the temperature and V-band veiling 
estimated by Stempels & Piskunov (2002), we calculate a photospheric luminosity L p h ot = 0.6 L & 
and R = 1.64 R & . Given these characteristics, the Siess et al. (2000) pre-main sequence stellar 
evolution tracks indicate that RU Lupi has a mass of 0.6 — 0.7 M Q and an age of 2-3 Myr. 

Although most previous studies have assumed that RU Lupi is a single star, Gahm, Petrov, 
& Stempels (2004) used the VLT to detect radial velocity changes of 2.5 km s _1 in photospheric 
absorption lines that are consistent with a period of either 3.7 or 5.6 days. They suggest that these 
changes are not due to spots, and may instead indicate the presence of a close stellar companion. 
RU Lupi is variable at ultraviolet and optical wavelengths. A photometric period of 3.7 d was 
also found by Hoffmeister (1958, 1965), but no period was confirmed in long-term monitoring (see 
Giovannelli et al. 1995). Hutchinson et al. (1989) and Giovannelli et al. (1990) found that RU Lupi 
is redder when it is fainter, which suggests that the brightening could be due to either an increase 
in the mass accretion rate, a decrease in circumstellar extinction, or cool starspots. 

Stempels & Piskunov (2002) used optical metallic absorption lines to measure a rotational 
velocity vsini = 9 ± 1 km s -1 . Rotational velocities of CTTSs can be at least vs'mi = 65 km s _1 
(Hartmann & Stauffer 1989; Basri k Batalha 1990; Johns-Krull & Valenti 2001). The measured 
rotational velocity suggests that either the rotational velocity is small or the rotation axis, and 
presumably the disk, are observed close to face-on. Ardila et al. (2002a) found that the equivalent 
width of the Mg II wind absorption in NUV spectra of RU Lupi is the largest of the 8 CTTSs in 
their sample, which we would expect if the wind is collimated and the disk is observed face-on. 
Blueshifted emission in profiles of the optical [O I], Na I, and Ca II H and K lines suggest that 
the disk obscures wind emission on the far side of the star, supporting the face-on viewing angle 
to the star (Stempels & Piskunov 2002). Based on the observations described above, we infer 
that the disk is probably observed close to face-on. If the 3.7 d period is due to rotation and 
not a close compation, the rotational broadening of absorption lines would imply an inclination 
of 24°. Spectro-astrometry of RU Lupi shows that Ha emission may be streaming SW and NE 
from the star (Takami et al. 2001). HH55, located 2' SW of the star, may be related to outflows 
from RU Lupi, which would imply that the disk is at least slightly inclined to our line of sight. 
The temperature and density of the base of the outflows from RU Lupi were also probed using 
HST/STIS observations of emission in Si III] 1892 A (G omez de Castro & Verdugo 2001). 

The ROSAT count rate of RU Lupi was ~ 9 x 10" 3 cts s 1 in a 15.8 ks PSPC observation 
(0.1-2.4 keV). We convert this to an X-ray flux of 8 x 10~ 14 erg cm" 2 s _1 , or L x = 5 x 1O~ 5 L 
using Neuhaiiser et al. (1995). 

RU Lupi is associated with a molecular cloud, mapped in 12 CO (J = 2—1) and 13 CO (J = 1—0) 
by Gahm et al. (1993) and in 12 CO (J = 2 - 1 by Tachihara et al. (2001). Gahm et al. (1993) 
found no direct evidence of any molecular outflows from RU Lupi, although they note a possible 
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connection with HH 55. Both studies found that the radial velocity of the molecular gas surrounding 
RU Lupi is about -0.4 km s _1 . 

3. OBSERVATIONS 
3.1. HST/STIS 

We observed RU Lupi with HST/STIS 1 as part of HST program 8157 (P.I. F. Walter). The 
spectra consist of five ~ 2500 s exposures with the E140M echelle grating, covering 1150-1710 
A, and a 120 s optical exposure with the G430L grating. We used the Of! 2 x / /06 aperture to 
minimize contamination by any off-source emission. The FUV observations use a 1024x1024 pixel 
Csl MAMA detector, and the optical observations use a CCD. The pixel size on the MAMA detector 
is 0'/036 in the dispersion direction (which corresponds to about 0.0142 A or 3.18 km s _1 at 1330 
A) and / /029 in the cross-dispersion direction. Details of these observations are listed in Table 1. 

The STIS FUV spectra were obtained in time-tag mode. Figure 1 shows that the count 
rate is similar at the beginning of each orbit but steadily increases within each orbit by about 
40%. Figure 2 shows that the width of the point-spread function in the cross-dispersion direction, 
measured in regions of strong line emission, decreased during these observations. We conclude 
that the telescope focus improved during each orbit. Thermal variations of HST can change the 
telescope focus. Models of thermal focus variations 3 indicate that the focus was poor at the start 
of each orbit but improved during every orbit. The nominal FWHM of the point-spread function in 
the FUV-MAMA is about 2.5 pixels. Because the size of the aperture (0 / /2 x 0706) in the dispersion 
direction is smaller than the line-spread function, the count rate increases as the focus improves. 

Several other observations obtained using the E140M grating and the / /2 x 0706 aperture on 
STIS, including those of V471 Tau, a white dwarf with a K2V companion (O'Brien, Bond, & Sion 
2001), and several CTTSs also suffer from this problem. On the other hand, observations of the 
two subluminous B stars PG1144+615 and PG1219+534 record a constant count rate during each 
orbit. No other STIS observations were taken in time-tag mode with the 072 x (y/06 aperture. This 
effect may also be seen data taken with small apertures that are not in time-tag mode. 

V471 Tau was observed at three different epochs, twice with the 0'.'2 x / /06 aperture and once 
with the 072 x 072 aperture (see Table 2, O'Brien et al. 2001; Walter 2004). The white dwarf 
dominates FUV emission and has a 555 s period (Jensen et al. 1986), during which its FUV flux 
varies by a few percent (Walter 2004). V471 Tau is also a bright FUV source, making it excellent 
to use for calibration. Figure 3 shows the total flux from 1250-1600 A in 300 s intervals of four 



2 See Leitherer et al. (2001) for a description of STIS. 

2 See http://www-int.stsci.edU/instruments/observatory/focus/focus2.html#breathing for a description of model- 
ing the HST focal length variations. 
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orbits from the observations of V471 Tau. We find one exposure taken with the small aperture, 
STIS observation number 24906, during which the detected flux is constant, at 4. Of x 10~ 10 erg 
cm -2 s _1 . Certain exposures, such as observation numbers 97314 and 97316, show flux increases at 
a rate about half that found in the observations of RU Lupi. These observations have a peak flux 
of 3.92 x 10 -10 erg cm" 2 s _1 . The observations taken with the large aperture have a flux 1250-1600 
A of 4.46 x 10" 10 erg cm" 2 s _1 . 

We measure the flux and cross-dispersion point-spread function in each 300 s interval for the 
observations of RU Lupi and V471 Tau. The point-spread function depends on wavelength, so 
we measure the FWHM of emission in the cross-dispersion direction at wavelengths with strong 
emission lines (O I, C II, Si IV, and C IV) in the spectra of RU Lupi, and in the same regions in 
the spectra of V471 Tau. Figure 4 shows the percent of the total flux between 1230-1650 A as a 
function of point-spread function for V471 Tau and RU Lupi. We calibrate the flux of RU Lupi by 
calculating a best-fit line to the data from V471 Tau, and then minimizing square of the difference 
between the data from RU Lupi and the best-fit line calculated from the V471 Tau observations. 

For RU Lupi, we calculate a total flux of 6.02 x 10 erg cm s" 1 between 1230-1650 A 
and 8.43 x 10" 13 erg cm -2 s" 1 between 1545-1555 A (a region dominated by C IV emission). 
We estimate that this flux calibration has an absolute uncertainty of about 15%. Inspection of 
the V471 Tau observations indicate that the changes in flux are constant with wavelength, such 
that the relative flux calibration across the wavelength range is uncertain to < 3% within a single 
exposure and between different exposures. The consistency of the relative count rate in strong lines 
during our observation of RU Lupi indicates that the relative flux is uncertain by < 10%. 

The Doppler correction in our data was in error by a factor of 1.6, such that the spectrum 
shifts by 3.8 pixels (12 km s" 1 ) during a single orbit. We corrected for this error by subdividing 
the spectra into 300 s intervals, and subsequently shifting and coadding the spectra. This error is 
prevalent in STIS echelle spectra in the FUV, but is now corrected in newer versions of the calstis 
pipeline. 4 The spectral resolution is about R = 45, 000. 

The wavelengths were initially calibrated using short exposures of the on-board Pt/Cr-Ne 
lamp obtained at the end of each observation of RU Lupi. We determined the absolute wavelength 
calibration independently by measuring the wavelength of the geocoronal Lya emission line. We 
then shift the line to 1215.58 A, where we expect the line to be at the time of our observation. 
We find that our wavelength solution corresponds well to that based on the calibration lamp. The 
absolute and relative wavelength solution is accurate to at least 3 km s _1 . 



4 see http: / /www. stsci.edu /hst /stis /calibration /pipe_soft_hist /update2 15c.html 
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3.2. FUSE 

We obtained a 25 ks spectrum of RU Lupi with the FUSE satellite in 21 separate exposures, 
each covering 912-1187 A, as part of FUSE program A109 (P.I. F. Walter). These observations 
occurred at night to minimize contamination by airglow emission. The FUSE has four co-aligned 
telescopes. The optical paths are coated with SiC or Al:LiF to maximize efficiency between 900- 
1000 A or 1000-1187 A, respectively. One observation produces eight overlapping spectra, each of 
which covers about 90 A. A detailed description of the FUSE mission can be found in Moos et al. 
(2000) and Sahnow et al. (2000). Since we used the LWRS aperture (30" x 30"), the FUSE spectra 
would include emission extended beyond the point source of RU Lupi, if present. 

The exposures were coadded and reduced with the calFUSE v2.4 software package (Kruk, 
Murphy, & Young 2001). The 8 separate spectra were then coadded using hrs-inerge in IDL. No 
significant differences were found in overlapping wavelength ranges of these segments. Since the 
LiF channels are more sensitive than the SiC channels, the S/N in the FUSE spectrum is higher at 
A > 1000 A. The data were binned by 4 pixels, which is smaller than the approximate instrument 
resolution of 9 pixels, or R ps 15,000. The flux scale is accurate to about ±5%, including a small 
degradation in some of the channels over time. 

We calibrate wavelengths in our FUSE spectrum by comparing the measured wavelengths of 
H2 emission lines (see §5.2) in the FUSE spectrum to those in the STIS spectrum. The H2 lines in 
the STIS spectrum can be characterized by a bright narrow component and a fainter and broader 
blueshifted component. When convolved to the lower spectral resolution of FUSE, the STIS H2 lines 
appear Gaussian with a velocity shift of -15.5 km s _1 relative to the photospheric radial velocity 
of RU Lupi. We offset the pipeline-calibrated wavelengths by +5 km s _1 , so that the average H2 
emission line in the FUSE spectrum has a velocity shift of -15.5 km s _1 . This velocity shift also 
places the peak of the N I 1135 A emission line at the radial velocity of the star, which is similar 
to the velocity shift of emission lines from neutral atoms in the STIS spectrum (see §5.2). The 
shift of 5 km s _1 is less than the 15 km s _1 accuracy of the standard calFUSE pipeline wavelength 
calibration (Kruk, Murphy, & Young 2001). We estimate an uncertainty of about 5 km s _1 in our 
wavelength calibration at A > 1100 A, and a larger uncertainty at A < 1100 A. 

The STIS spectrum overlaps the FUSE spectrum between 1160-1185 A, although the sensi- 
tivity of STIS is poor at these short wavelengths. Figure 5 shows that the C III line at 1175 A is 
much stronger in the STIS spectrum than in the FUSE spectrum, so fluxes from the two spectra 
are not directly comparable. 
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4. THE ISM IN OUR LINE OF SIGHT TO RU LUPI 

4.1. Interstellar Absorption Lines 

Narrow absorption lines, produced in the interstellar medium (ISM), are detected against 
several stellar emission lines (see Table 3 and Figure 6). In §5.2, we find that most low-ionization 
emission lines have Gaussian profiles centered about the radial velocity of RU Lupi. We therefore 
assume that the intrinsic emission profile is Gaussian to measure the velocity and equivalent width 
of the absorption. 

Interstellar absorption is centered between —8 and —19 km s _1 in the heliocentric frame, 
depending on the line. RU Lupi is located at / = 338°, b = +11°89, where models of the local 
interstellar medium produced by Redfield & Linsky (2000) predict that very little absorption by 
the LIC is expected and some absorption by the G cloud is possible. These models also predict that 
any absorption in the local interstellar cloud and the G cloud should be centered at -23 km s" 1 and 
-26 km s _1 , respectively. The Lupus molecular cloud has a velocity of about -0.4 km s _1 in our line 
of sight toward RU Lupi (Gahm et al. 1993; Tachihara et al. 2001). Tachihara et al. (2001) estimate 
that the H 2 column density of the Lupus cloud in the direction of RU Lupi is log N(R 2 ) = 21.0, 
which is much higher than the column density 19.0 < logiV(H2) < 19.5 in our line of sight to RU 
Lupi, suggesting that the star is located in front of most of the cloud. The data show absorption 
components redshifted relative these clouds, indicating that other interstellar clouds are located 
along the RU Lupi line of sight. 

Ardila et al. (2002b) found that the equivalent widths of interstellar Mg II lines towards 
CTTSs are larger than typical values measured through the ISM toward nearby sources. We also 
detect large equivalent widths in many narrow absorption lines in the FUV, including a detection 
of interstellar absorption in the Si III 1206.5 A line. Interstellar Si III absorption has only been 
detected in one other local sightline, towards e CMa (Gry et al. 1995; Wood et al. 2002b). The large 
Mg II equivalent widths and detection of Si III absorption may be produced by warm or shocked gas 
in the Lupus molecular cloud. Any additional absorption components due to the Lupus molecular 
cloud, a remnant envelope, or a slow wind should be centered within a few km s _1 of the interstellar 
absorption. 

Lya is detected in emission only longward of 1217.5 A, indicating that the sum of interstellar 
and circumstellar H I absorption is optically thick shortward of 1217.5 A. We do not detect any 
blueshifted emission from Lya because of strong H I absorption in the wind. ISM absorption should 
dominate absorption on the red side of the line. We model the observed Lya profile with a single 
Gaussian emission profile, centered at the radial velocity of the star, and a Voigt absorption profile. 
The height of the Gaussian is determined from the best fit to the red side of the Lya line for a 
given iV(H I). We assume that the absorption is centered at -15 km s _1 in the heliocentric frame, 
roughly the velocity at which we detect interstellar absorption in lines such as O I and C II (Table 
3). We assume a Doppler broadening parameter b ~ 11 km s _1 , assuming a temperature of 7000 
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K, typical of the local ISM (Redfield k, Linsky 2000). At the high column densities applicable to 
the RU Lupi sightline, the column density is insensitive to b. Figure 7 compares the profiles that 
best fit the observed emission for a range of neutral hydrogen column density 7V(H I). We find that 
log iV(H I)= 20.0 ± 0.15 (cm" 2 units), where the error includes an estimate of the uncertainty in 
the central velocity of the ISM absorption. 

The low S /N in the continuum prevents us from detecting most of the H2 absorption lines that 
should occur in the FUSE bandpass. However, the 6-0 P(3) line at 1031.19 A is detected because 
it occurs against strong O VI emission (Fig. 8). This H2 line is shifted by about —23 ± 5 km s _1 
relative to the stellar radial velocity. The O VI 1037 A line is weak and narrow because most of 
the line is absorbed by H2 and C II (see Fig. 8). The detection of the one H2 line and only weak 
O VI 1037 A emission imply that 19.0 < logiV(H 2 ) < 19.5 for temperatures 100 K< T(H 2 ) < 300 
K, which is typical of interstellar molecular hydrogen (Rachford et al. 2002). We therefore measure 
a total hydrogen column density of log iV(H)=log[iV(H I) + 2 x 7V(H 2 )] = 20.1 ± 0.2. Stempels 
& Piskunov (2002) measured optically thin absorption in the Na I line at 5890 A, and found 
logiV(Na l)= 11.0. They used interstellar abundances from Welty & Hobbs (2001) to estimate 
logiV(H)~ 19.8, which is approximately what we measure. 



4.2. Dust Extinction 

Strong optical veiling due to accretion prevents a reliable estimate of the extinction using 
photospheric emission from RU Lupi. Instead, we use the interstellar absorption discussed above to 
constrain the extinction. Assuming a typical interstellar gas/dust ratio and dust grains described 
by Ry = 3.09 (Cardelli, Clayton, & Mathis 1989), we convert the hydrogen column density of 
log 7V(H)= 20.1 ± 0.2 to an extinction Ay ~ 0.07 mag using the relation described by Bohlin et al. 
(1978). This low extinction implies that RU Lupi is near the front edge of the Lupus molecular 
cloud. We adopt Ay = 0.07 mag for the remainder of this paper. We confirm a low extinction to 
RU Lupi using two other probes in our data. 

The NUV and blue continuum of RU Lupi is produced by accretion (e.g. Hartigan et al. 
1991). We use the shape of the accretion continuum to constrain the extinction by modelling the 
continuum as emission from an isothermal slab, described in §5.1. Model continua that best fit the 
observed UV-optical (1650 A< A < 5700 A) continuum are calculated for each assumed value of the 
extinction. Figure 9 shows the estimated accretion continuum, calculated independently for three 
assumed extinctions compared with the extinction-corrected NUV spectrum of RU Lupi. When 
we use Ay > 0.4 mag, our models of the accretion continuum underestimate the flux at 2175 A 
because the extinction bump at 2175 A increases the assumed intrinsic flux at that wavelength. 

If instead most of the dust grains in our line of sight to RU Lupi are in the Lupus molecular 
cloud, then the dust may differ from typical interstellar dust. The grain size should not differ 
significantly from grains typical of the local ISM. Although sightlines through the densest regions 
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of molecular clouds can have Ry = 5.5, Kenyon, Dobrzycka, & Hartmann (1994) and Whittet et al. 
(2001) find that the sightlines through the Taurus and Ophiucus molecular clouds with Ay < 3.0 
do not vary significantly from the average interstellar value of Ry = 3.09 (Cardelli et al. 1989). 
However, the extinction law towards the peculiar B star HD29647 {Ay = 3.6 mag), which lies 
behind Taurus, does not show the strong 2175 A bump (Whittet et al. 2004) that is typically 
associated with small carbonaceous dust grains (e.g. Draine k, Malhotra 1993). The grains in the 
line of sight to HD29647 have an Ry = 3.63, which suggests that they are large relative to typical 
interstellar grains, described by Ry = 3.09, (Cardelli et al. 1989). Calvet et al. (2004) applied 
the HD29647 extinction law to analyze ultraviolet spectra of intermediate-mass T Tauri stars in 
Taurus. If the extinction law of HD29647 applies to low-extinction targets such as RU Lupi, then 
the shape of the accretion continuum is not a good probe of extinction because of the absence of 
the 2175 A bump. 

We also place an upper limit on the extinction by comparing the relative strength of fluorescent 
H2 emission lines (see §5.2), which occur between 1100-1650 A, to the same emission lines from 
TW Hya. Extinction is negligible to TW Hya, and both stars are observed nearly face-on, so any 
geometrical effects on the H2 fluxes should be similar. In §5.2, we identify 90 unblended H2 lines in 
the STIS spectrum of RU Lupi from 15 different upper levels. In the STIS spectrum of TW Hya, 
Herczeg et al. (2002) detected 146 H2 lines from 19 different upper levels. By comparing the relative 
flux in lines from the same upper level for the two stars, we place an upper limit of Ay < 0.6 mag 
toward RU Lupi for either a typical interstellar extinction law or the interstellar extinction law of 
HD29647. 

Our estimate of a low extinction towards RU Lupi agrees with most previous estimates. Gahm 
et al. (1979) estimated 0.3 < Ay < 1.0 mag using the IUE spectrum of RU Lupi. Giovannelli et 
al. (1995) estimated an extinction of Ay = 0.3 mag based on the 2200 A reddening bump in IUE 
spectra, and their result agrees with the extinction Ay = 0.4 mag estimated by Brooks et al. (2001) 
and Lamzin et al. (1996). Our extinction estimate is much lower the extinction of Ay = 1.28 mag 
calculated by Hughes et al. (1994) based on optical colors, which are contaminated by emission due 
to accretion. 

The molecular fraction of gas towards RU Lupi is between 17-40%. Most sightlines with a 
molecular fraction > 0.01 have Ay > 2.4 mag (Savage et al. 1977; Rachford et al. 2002). However, 
this relationship may not apply to environments around young stars. A high molecular fraction 
and low extinction were also found towards the Herbig Ae star HD104237 (Grady et al. 2004). 
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5. ACCRETION 

5.1. Mass Accretion Rate 

We use the accretion luminosity, estimated from the UV excess and the Balmer jump at 3700 
A, to approximate the mass accretion rate of RU Lupi (Valenti et al. 1993; Gullbring et al. 1998). 
Our HST observations consist of nearly simultaneous optical and FUV spectra. We obtain NUV 
spectra of RU Lupi using a combination of archival observations from HST /STIS, ffi>T/GHRS, 
and IUE (Table 4). The NUV spectra are normalized by factors listed in Table 4 to match the 
optical and FUV continua flux, ignoring regions of strong emission lines. We scale the optical STIS 
spectrum of V836 Tau, a K7 WTTS, to use as a template for the photospheric emission from RU 
Lupi. We estimate that Ay = 1.52 mag towards V836 Tau, the average extinction estimate from 
Kenyon & Hartmann (1995) and Whittet et al. (2001). 

The methodology used to calculate the mass accretion rate is described by Valenti, Basri, 
& Johns (1993) and Herczeg et al. (2004). We model the spectrum of the accretion shock by 
approximating the emission from an isothermal, pure-hydrogen slab (see Fig. 10), which fit the 
data well. Other studies have also calculated the accretion continuum by modelling the accretion 
shock (Calvet & Gullbring 1998). The emission spectrum that best fits the observed spectra is 
produced by a slab with temperature T = 9600 K, electron density n e = 3 x 10 14 cm" 3 , slab 
surface area parameter (fraction of stellar surface area covered by accretion) 5 = 0.017, thickness 
/ = 2 x 10 7 cm, and a scaling factor for the template photospheric spectrum of 0.39. The slab has 
an optical depth of 0.64 at 5000 A, which is less than the average optical depth in slab models of 
blue emission from other CTTSs (Basri & Bertout 1989; Hartigan et al. 1991). This scaling factor 
corresponds to an optical veiling in the continuum of a factor of 4-8 between 4000-4300 A, and 
about 2 at longer wavelengths. We calculate an accretion luminosity = 0.35 L . We then 

estimate a mass accretion rate of 

M . = LggAg R 0.8M e 

GM t V ' 140 2 1.65^0 M y ' 

This mass accretion rate is large, but significantly smaller than previous rates that were estimated 
assuming Ay = 1.3 (e.g., Hughes et al. 1994; Ardila et al. 2002a). Although mass accretion rates 
can be as high as 10~ 6 M yr _1 for CTTSs, the mass accretion rate of 5 x 10" 8 M Q yr _1 is high 
for CTTSs with ages of 2-3 Myr (Hartmann et al. 1998). 

The model accretion continuum accounts for the observed optical and NUV continuum. A 
slight FUV excess may be detected, although the model accretion continuum in the FUV is sensitive 
to uncertainties in the temperature of the accreting gas. Figure 11 shows that the FUV continuum 
from RU Lupi, if present, is much weaker than that detected towards TW Hya and several other 
CTTSs (Herczeg et al. 2004; Bergin et al. 2004). The FUV continuum from CTTSs is most likely 
produced by H2 fluorescence generated by highly energetic electrons colliding with and electronically 
exciting H2 molecules, and may be related to an inner dust hole in the disk (Liu et al. 2002; Bergin 
et al. 2004). A bump at 1600 A, diagnostic of this H2 emission, is weakly present. 
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5.2. FUV Emission Lines 

The total observed emission between 1170-1700 A from RU Lupi is 7.16 x 10~ 12 erg cm 
s _1 , most of which occurs in emission lines. The flux in hot FUV emission lines from CTTSs 
increases with the mass accretion rate (Johns-Krull et al. 2000; Calvet et al. 2004). Many Fe II 
and fluorescent H2 lines are pumped by other emission lines, such as Lya, so their strength should 
also scale with the accretion rate. Johns-Krull et al. (2000) found that H2 flux scales with the C IV 
flux. Other lines, such as the O I 1305 A triplet and the C II 1335 A doublet, may be produced 
both in accreting gas and in the wind. Thus, we can analyze certain emission lines to probe the 
accretion flow. 

For emission lines that are not corrupted by wind absorption, we fit a single or multiple 
Gaussians to the emission line profile to measure the line center, FWHM, and flux in the line. 
Tables 5 and 6 list parameters for identified atomic or unidentified emission lines detected in the 
STIS and FUSE spectra, respectively. For lines that have large wind absorption components, we 
list only the continuum-subtracted flux detected in the line. Fluxes are not directly comparable 
between the STIS spectrum and the FUSE spectrum owing to temporal variability in the emission 
(see Fig. 5). 

Strong resonance lines, such as the C IV doublet at 1549 A, are commonly detected in emission 
from late-type stars and are also strong in the FUV spectrum of RU Lupi. We identify weak emission 
lines using the atomic spectral line database of Kurucz & Bell (1995). The identifications marked 
by a '?' in Tables 5-6 are uncertain. Since many lines, including fluorescent H2 emission, Fe II, 
and Cr II, are pumped by Lya, our identification technique relies on searching for lines that could 
be pumped by Lya or other strong emission lines. We list the pumped lines and the pumping 
transition in Table 7. Many of these lines were identified by Johansson & Jordan (1984), Carpenter 
et al. (1988), and Hartman & Johansson (2000). Other pumped lines include CI I 1351.6 A, which 
is pumped by C II 1336 A, and S I 1472 A, which is pumped by O I 1302 A and S I 1295 A. The 
O I triplet at 1304 A can be pumped by Ly/? 1026 A (Carpenter et al. 1988). Although airglow 
Ly/3 prevents us from detecting the center of Ly/3 in our FUSE spectrum, we do not detect Ly/3 
emission at 100 km s _1 . Nevertheless, the O I lines are very strong. Moreover, other lines such as 
C II that require similar excitation are also very strong. We speculate that the O I emission in the 
RU Lupi spectrum is most likely produced by collisional excitation or recombination in the hot gas 
rather than by fluorescence. 

Many narrow (FWHM~ 17.5 km s _1 ) emission lines detected in the STIS spectrum are identi- 
fied as H2 lines (see Table 8) based on the H2 transition database calculated by Abgrall et al. (1993) 
and the linelist from #ST/STIS observations of TW Hya (Herczeg et al. 2002). Most of these lines 
are Lyman-band (B-X) H2 lines pumped by Lya. We also detect one Werner-band (C-X) H2 line 
at 1208.8 A that is pumped by O VI emission. These H2 lines are found in UV spectra of CTTSs 
(Brown et al. 1981; Valenti et al. 2000; Ardila et al. 2002a), red giants (Ayres, Brown, & Harper 
2003a), and the wind of Mira B (Wood et al. 2002a). We coadded the strongest H 2 lines in the STIS 
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spectrum to analyze the line profiles. Figure 12 shows that the coadded H2 lines are well fit by a 
strong Gaussian, centered at -12 km s -1 , and a weaker Gaussian, centered at -30 km s _1 , relative 
to the radial velocity of RU Lupi. We also detect additional flux in these lines from v = —120 to 
—70 km s _1 . The flux ratio of these two Gaussians is 1.3 and the ratio of the widths is 0.33. To 
measure fluxes, we fit each individual H 2 line with two Gaussians, with the velocity separation, 
relative width, and relative flux fixed to the numbers obtained from the fit to the coadded profile 
above. We also detect unresolved Lyman and Werner-band H2 lines with FUSE. We fit a single 
Gaussian profile to these lines, and list their fluxes in Table 9. If present, any extended H2 emission 
would be detected by FUSE because it has a large (30 / . / x 30'/) aperture. Figure 13 compares the 
coadded spectral line profiles of the H2 lines in the STIS and FUSE spectra, and the geocoronal 
Ly/3 line in the FUSE spectrum. The FUSE H2 lines are not significantly broader than the STIS H2 
lines, convolved to the lower resolution of FUSE. If the H2 emission were extended beyond a point 
source in the dispersion direction, we would expect the spectral line profile would be much wider. 
This comparison therefore suggests that the H2 emission in the FUSE spectrum is not significantly 
extended beyond a point source. The H2 lines will be discussed in detail in a future paper. 

Some cooler emission lines typically found in stellar chromospheres, such as the O I 1304 A 
triplet and the C II 1335 A doublet, are partially or completely absorbed by wind features listed in 
Table 10 (see also Fig. 16). Figure 14 shows that emission in these cool lines are centered at the 
radial velocity of RU Lupi. 

Emission in C IV, Si IV, N V, and possibly O VI is most likely produced by 10 5 K gas heated 
by the accretion shock. Figure 15 shows that these lines are strong and broad. The STIS lines 
have centroids near the radial velocity of the star. The FUSE lines of O VI 1032 A and C III 977 A 
lines are blueshifted by ~ 190 and ~ 100 km s _1 , respectively. Low S/N prevents any detection of 
wind in the C III 977 A line, but if present we may be underestimating the blueshift. Absorption 
in H2 and other lines may mask some redshifted emission in the O VI and C III lines. The velocity 
shift of these lines may be time dependent. Alternatively, the emission in C III and O VI resonance 
lines during the FUSE observation might be produced by shocked emission due to the outflow, as 
hot gas can be seen in some HH objects (Raymond et al. 1997). The lack of C III 1175 A emission 
suggests that any emission produced by accreting gas is faint during this observation. FUSE would 
then only detect this spatially extended blueshifted emission. The STIS observations would most 
likely not be sensitive to this extended emission because of its small aperture. The time dependence 
in the strength of C III emission is discussed in §7. 

The C IV 1548 A line shows broad blueshifted emission, extending to -300 km s _1 , that is 
not present in the other emission lines. In contrast, the Si IV 1394 A line has stronger redshifted 
emission than the Si IV 1403 A line. The difference in behavior between the C IV doublet and the 
Si IV doublet is puzzling. The N V 1243 A line is absorbed by N I in the wind (Fig 16), while 
interstellar H2 absorption corrupts the O VI line profiles (Fig. 8). Anomalously weak N V 1243 A 
emission was also detected in the FUV spectrum of T Tau (Walter et al. 2003). 
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Optical depth effects likely produce these complex emission features. The mass accretion rate 
is large, so we infer that the accretion column may be optically thick. The C IV 1548 A and 1550 A 
lines from TW Hya have a similar shape to each other (Herczeg et al. 2002), possibly because the 
mass accretion rate of TW Hya is 25 times smaller than that of RU Lupi, so the accretion shock 
may be less optically thick. We also note that, unlike TW Hya, Si II, Si III, and Si IV lines from RU 
Lupi are all bright. Understanding these profiles will require detailed modelling of the accreting 
gas. 

Most of the observed Lya emission is absorbed by H I in the stellar outflows, in the ISM, 
and possibly the accretion shock itself. In Herczeg et al. (2004), we found that the intrinsic Lya 
emission could be reconstructed from the H2 lines, and may contribute 85% of the total FUV 
emission from CTTSs. The strength of the H2 and Fe II lines that are pumped by Lya imply that 
the Lya emission from RU Lupi is similarly strong. We will estimate the Lya emission from RU 
Lupi in a future paper. 

Coronal Fe lines have been detected in UV spectra of the Sun (Sandlin, Brueckner, & Tousey 
1977) and cool stars (Ayres et al. 2003b; Redfield et al. 2003). Lamzin et al. (2001) predict that 
emission in the [Fe XII] 1349.382 A and Fe XI 1467.420 A lines may be produced at the accretion 
shock and could be detected in UV spectra of CTTSs. Based on models of the accretion shock by 
Lamzin et al. (2001), the lines should be centered at v s hock/^- The FWHM of these emission lines 
should be dominated by thermal broadening. We do not detect any coronal Fe lines in the RU Lupi 
spectra, and place unreddened upper limits of ~ 2 x 10 -15 erg cm" 2 s _1 that are about an order 
of magnitude larger than the line fluxes predicted by Lamzin et al. (2001). 

6. WIND ABSORPTION 

Outflows from RU Lupi are detected by the strong blueshifted absorption seen in P-Cygni 
profiles and against the continuum (Fig. 6). The absorption is typically seen from low energy levels 
of neutral or singly-ionized species. 

The cool emission lines that are not corrupted by wind absorption are symmetric about the 
radial velocity of the star (see §5.2 and Figure 14). We therefore estimate the intrinsic profile of 
lines with wind absorption by reflecting the red emission about the radial velocity of the star. We 
then measure the equivalent width and central and maximum velocities of the absorption feature 
(see Table 10). Lines such as O I, C II, and Fe II have P-Cygni profiles, with absorption ranging 
from about —70 to —270 km s -1 . For certain lines of Fe II and N I, the wind absorption is detected 
against the continuum. The wind absorption feature seen in Si III is narrower than the wind 
absorption features detected in other lines from the ground state, which may be due either to a 
small column density of Si III or the absence of a velocity gradient for the hot gas in our line of 
sight. The latter possibility could occur if most of the Si III occurs close to the star. Most of the 
wind absorption is detected in transitions with low-energy lower levels. However, we also detect 
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wind absorption in the N I lines at 1242, 1492, and 1494 A from lower levels excited to 2 eV (Fig. 
16). 

In certain cases, the wind absorption is strong enough that many nearby lines are barely 
detected. Figure 17 shows the O I triplet at 1304 A and the C II doublet at 1335 A. The wind in 
the C II 1335.7 A line obscures emission in C II 1334.5 A. Similarly, the O I 1305 A line is weak 
due to wind in that line and a wind in O I 1306 A. 

7. VARIABILITY 

Figure 5 shows a striking decrease in flux in the C III 1175 A multiplet between the STIS and 
the FUSE observations. In the FUSE spectrum, most of the emission detected near 1175 A is H 2 
and not C III. The O VI 1035 A doublet and C III 977 A emission are weak and blueshifted. If 
the drop in flux seen in C III 1175 A is related to a decrease in the mass accretion rate, then we 
may also not be detecting any emission produced by the accretion shock in the VI 1035 A doublet 
or the C III 977 A line. Instead, with its large aperture, FUSE may have detected emission in 
these lines that is produced by shocks in the outflow. The fluorescent H2 emission in the FUSE 
spectrum, which is highly sensitive to the Lya emission, is comparable to the flux we expect based 
on the STIS spectrum. In §5.2, we found that the H2 emission in the FUSE bandpass originates 
predominantly on-source. Because the H2 lines serve as a proxy for the strength of Lya emission, 
we infer that the Lya emission has a similar flux in the two observations. This result is surprising 
because we expect that both C III and Lya emission are produced by accreting gas. We speculate 
that a smaller mass loss rate significantly reduces the optical depth of the C III line, while the Lya 
line remains optically thick. The FUV continuum where the STIS and FUSE spectra overlap is 
marginally detected at ~ 2 x 10 erg cm" - 2 s- 1 A" 1 . This emission, which is likely produced by 
H2 (Bergin et al. 2004), did not change by more than ~ 50%. 

We also compare the observed STIS spectrum with spectra obtained using the G160M grating 
on HST /GHRS, which cover narrow (35 A) spectral regions centered on the C IV and Si IV lines 
(Ardila et al. 2002a). In the GHRS spectra, the FUV emission lines have similar shapes but are 
four times weaker than the fluxes detected with STIS (Fig. 18). The continuum in the GHRS 
observation also appears weaker than that detected in the STIS observation. However, the H2 
emission line fluxes are similar in the STIS and the GHRS spectra, which suggests that the Lya 
emission did not change significantly between the two observations. The GHRS observation could 
also include off-source H2 emission, as is detected towards T Tau (Walter et al. 2003). Our STIS 
observations must have fortuitously occurred when the accretion rate of RU Lupi was higher than 
usual. 

We find no significant variability during the STIS spectrum. Although this analysis is compli- 
cated because the count rate depends on the changing focus of HST (see §3.1), we place a limit that 
the flux of RU Lupi varied by no more than 10% in intervals longer than 100 s and 4% in inter- 
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vals longer than 300 s. Some minor differences appear in the profiles of hot lines during the STIS 
observations. The S/N of the FUSE observations prevented any significant analysis of variability 
over short timescales. 

8. CONCLUSIONS 

We present and analyze FUV spectra of RU Lupi, a 2-3 Myr old CTTS with a high mass 
accretion rate. Based on these spectra, we find the following results: 

1. The HST telescope focus improved during each orbit of the RU Lupi observation, which 
increased the detected count rate because the FWHM of the line-spread function was larger than 
the aperture size in the dispersion direction. We account for this instrumental problem by using 
observations of V471 Tau. This problem appears in several other time-tag STIS spectra that used 
a small aperture. 

2. We estimate a total hydrogen column density log iV(H)= 20.1 ±0.2 in our line of sight to RU 
Lupi using Lya and H2 absorption, and convert the column density to an extinction of Ay ~ 0.07 
mag. 

3. The mass accretion rate of RU Lupi is about (5 ± 2) x 10~ 8 M Q yr _1 . 

4. Emission is detected in 174 lines, 103 of which are Lyman and Werner-band H2 lines pumped 
by Lya or O VI. Some Fe II and Cr II lines are also pumped by Ly<x The hot lines of C IV, Si IV, 
and N V are bright and broad, but differences in the line profiles are puzzling. In particular, the 
C IV 1548 A line has strong redshifted emission that is not seen in the C IV 1550 A line, but the 
Si IV 1394 A line has strong blueshifted emission relative to the Si IV 1403 A line. Emission in the 
O VI 1035 A doubet and C III 977 A is blueshifted, and at the time of our FUSE observation may 
be dominated by emission produced by shocks in the outflow. The observed Lya emission only 
accounts for about 1% of the observed FUV emission, likely because of strong absorption by H I in 
the ISM and stellar wind. 

5. Wind absorption is detected in 30 different lines, and generally ranges from -70 to -270 
km s _1 . Atoms with low first-ionization potentials are ionized in the wind, resulting in absorption 
in lines of Si II, Si III, Fe II, Al II, and C II. Wind absorption is also detected in lines of N I and 
O I. Three N I wind absorption lines originate from levels excited 2.38 eV above the ground state. 

6. The FUV emission from RU Lupi varies significantly over months-long timescales. Our 
STIS observation of RU Lupi occurred during a period of strong emission, possibly due to a high 
mass accretion rate. The C III emission is much weaker in the FUSE spectrum than in the STIS 
spectrum, although the H2 emission, which is a proxy for Lya emission, does not vary significantly 
between these two observations. 
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Fig. 1. — In our STIS observations of RU Lupi, the count rate within each orbit increased with time, 
but then returned to a low level at the start of the next orbit. This increase of flux results because 
the FWHM of the line-spread function is wider than the aperture in the dispersion direction but 
narrows during each orbit, due to thermal changes in the telescope. 
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Fig. 2. — The FWHM of C IV and Si IV emission in the cross-dispersion direction, measured in 
300 s intervals, improves during each orbit because of thermal focus variations. 
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Fig. 3.— The flux detected between 1250-1600 A in four of 16 orbits of STIS observations of V471 
Tau. A gradual increase of detected flux with time is present in two of the orbits shown here, and 
is similar to the increase seen in our spectra of RU Lupi. 
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Fig. 4. — The flux between 1230-1650 A from RU Lupi (diamonds) and V471 Tau (asterisks) 
depends on the telescope focus. The fraction of the total flux is calculated by comparing the 
measured flux in a 300 s interval to the estimated total flux, for the wavelength region between 
1230-1650 A. The image size in the dispersion direction is plotted here as the average of the 
FWHM of O I (1292-1313 A), C II (1333-1337 A), Si IV (1392-1405 A), and C IV (1547-1552 
A). The dotted line is the line that best fits the data from the calibration source V471 Tau. We 
then estimate the total flux from RU Lupi by fitting the RU Lupi data to this line using the least 
squares minimization method. The point-spread function is broader than the small aperture, so 
that better focus improves the count rate on the detector increases. 
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Fig. 5. — C III emission was strong during the HST/STIS observation (dotted line), but was much 
weaker during the FUSE observation (solid line), which occurred about a year later. This striking 
difference likely resulted because of a drop in the mass accretion rate between the two observations. 
The wavelengths of fluorescent H2 emission lines are indicated by the dashed vertical lines. 
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Fig. 6. — P-Cygni profiles are detected for many lines, four of which are shown here. The wind 
absorptions in most lines range from ~ —70 — 270 km s _1 , although the absorption in the Si III 
line is much narrower than wind absorption from ground state transitions of other species because 
the line is less optically thick. Narrow absorption near line center (except for the excited O I 1306 
A line) is produced by the interstellar medium and possibly also by material local to RU Lupi. 
Interstellar Si III absorption is not detected locally (Wood et al. 2002b). The dashed vertical line 
marks the average velocity of interstellar absorption detected towards RU Lupi. 
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Fig. 7. — The observed red wing of the Lya line (solid line) is modeled as a Gaussian profile, 
centered at the heliocentric velocity of Lya, and a Voigt absorption profile for interstellar H I, 
centered at the interstellar velocity measured from other absorption lines (see Table 5). No Lya 
emission is detected on the blue side of line center. We show here the fit for three different assumed 
column densities of H I. We determine that log iV(H 20.0 ± 0.15. 
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Fig. 8. — The FUSE spectrum of RU Lupi (solid line) is compared with the transmission percentage 
(dashed line) through gas with T = 300 K and log iV(H2)=19.0. Both lines in the O VI doublet are 
blueshifted from line center (dotted vertical lines). The O VI 1032 A profile shows strong emission, 
with absorption in the H2 line at 1031.19 A. The O VI 1037 A line is detected, but most of the line 
is absorbed by C II (dashed vertical lines) and H2 in the line of sight to the star. 
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Fig. 9. — We fit models (dashed lines) of the accretion continuum to the observed UV-optical 
spectrum of RU Lupi, assuming a range of extinctions. When we assume Ay > 0.4 mag, we cannot 
reconcile the increased emission at 2175 A with the emission at shorter wavelengths. Low extinction 
is also consistent with the observed N(R I) column density in our line of sight to the star (see §4.1). 
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Fig. 10. — We model the observed accretion continuum (solid line) with emission from a pure 
hydrogen slab (dashed line). The accretion continuum is obtained from data described in Table 3, 
assuming Ay = 0.07. To extract the accretion continuum shown here, we subtracted the scaled 
spectrum of the WTTS V819 Tau, a template for the photospheric emission of RU Lupi, from the 
observed spectrum of RU Lupi. 
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Fig. 11. — Comparison of the observed UV emission from RU Lupi (bottom solid line) and TW Hya 
(top solid line, multiplied by 10) with spectra calculated from modelling the accretion continuum 
(dashed lines). The FUV continuum in excess of emission from the accretion column is much weaker 
from RU Lupi than from TW Hya. 
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Fig. 12. — The spectral profile of coadded H2 lines from the STIS spectrum (solid line) shows that 
the lines are blueshifted and asymmetric. We fit the profile with two Gaussian profiles (dotted). 
The summed profiles are shown by the dashed line, with parameters described in §5.2. Excess 
blueshifted emission at v = —120 to —70 km s _1 is also detected. These lines will be discussed in 
detail in a future paper. 
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Fig. 13. — We compare the spectral profile of coadded H2 lines observed with FUSE (solid line) 
with geocoronal Ly/3 emission (dotted line) and coadded H2 lines from the STIS spectrum (dashed 
line), convolved to the spectral resolution of FUSE. The geocoronal Ly/3 emission fills the aperture 
and in this plot is shifted in velocity to the center of the H2 lines. The coadded FUSE H2 lines 
have a similar FWHM as the coadded STIS H2 lines, which suggests that the H2 emission detected 
with FUSE is dominated by on-source emission. Any emission extended in the dispersion direction 
would broaden the line profile. 
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Fig. 14. — The weighted mean of the radial velocity of cool (neutral, singly, and doubly ionized) lines 
is -1.4 km s _1 (dotted line) relative to the radial velocity of RU Lupi. Lines that are contaminated 
by wind emission and lines that are only tentatively identified (marked by a ? in Table 5) are not 
used in this analysis. 



-34- 




Velocity (km s ! ) 



Fig. 15. — Emission lines produced in ~ 10 5 K gas are broad and strong in spectra of CTTSs. These 
lines from RU Lupi show complicated features. The C IV 1548 A line shows a strong blue bump 
that is not detected in any of the other lines, while the Si IV 1394 A line shows excess redshifted 
emission relative to the other lines. N V 1243 A emission is absorbed by N I in the wind (see Fig. 
16). The O VI and C III lines detected in the FUSE spectrum appear blueshifted and may be 
produced by shocks in the outflow. Because the FUSE observation occurred about one year after 
the STIS observation and the FUV emission from RU Lupi is variable, the O VI and C III lines 
shown above should not be directly compared with the C IV, Si IV, or N V lines. 
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Fig. 16. — Wind absorption is detected in N I lines from excited lower levels, with E=2.38 eV. The 
N I absorption near 1243 A attenuates the N V 1243 A emission line. 
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Fig. 17. — Wind absorption in the O I and C II lines. The wind absorbs a significant amount of the 
emission line and also emission from adjacent lines. O I 1305 A is weak because of wind absorption 
due to O I 1306 A. Similarly, C II 1334.5 A is barely detected because of absorption in C II 1335.7, 
and two H2 lines are weak because of absorption in C II 1334.5. Line wavelengths at the radial 
velocity of RU Lupi are indicated by the vertical dotted lines. The narrow absorption feature near 
line center is produced by the ISM. 
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Fig. 18. — The STIS spectrum of RU Lupi has a stronger FUV continuum and stronger emission 
lines than the previously observed GHRS spectrum (Ardila et al. 2002a). 
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Table 1: OBSERVATIONS OF RU LUPI 

Date Instrument Exposure (s) Bandpass (A) Aperture Spectral Res. (A/AA) 

7/12/00 a HST /STIS E140M 12530 1150-1700 / /2 x 0706 45000 

7/12/00 a HST /STIS G430L 120 2900-5700 52'/ x 071 800 

9/28/01 6 FUSE 24517 900-1187 30" x 30" 15000 

a JD 2451738 6 JD 2452180 

Table 2: HST/STIS OBSERVATIONS OF V471 TAU 

Dates Orbits Exposure Aperture ^"400 Program Obs. Numbers 6 P.I. 

3/13-3/23/98 c 6 14350 0"2 x 0"06 401 7735 24906-24919 F. Walter 

8/24/00-8/27/00^ 4 7848 072 x 0"2 4.46 8305 73220-73229 H. Bond 

l/25/02 e 4 10336 0"2 x 0706 4.01 9283 97310-97316 F. Walter 

a Peak flux from 1250-1600 A, in 10" 10 erg cm s STIS observation numbers 

C JD 2450885-2450896 d JD 2451781-2451784 e JD 2452299 



Table 3: Properties of Interstellar Absorption Lines" 



ID X obs (A) v (kms- 1 ) W\ (A) 



Si III 


1206.447 


-11.2 


0.089 ± 0.007 


O I 


1302.090 


-16.2 


0.186 ±0.008 


CII 


1335.662 


-8.3 


0.098 ± 0.005 


Sill 


1526.618 


-15.1 


0.172 ±0.009 


Fell 


1608.367 


-13.7 


0.107 ±0.010 


Al II 


1670.672 


-18.7 


0.173 ±0.008 



"Measured by assuming Gaussian emission lines. 



Table 4: SPECTRA OF THE ACCRETION CONTINUUM 



Bandpass 


Instrument 


Scaling Factor" 


S/N fe 


Date 


1185-1700 


HST/STIS E140M 


1.0 


20 


07/12/00 


1700-1839 


IUE SWP-LO 


0.75 


60 


many c 


1839-1930 


HST/STIS G230M 


1.1 


35 


08/24/92 


1930-2302 


IUE SWP-LO 


0.75 


14 


many c 


2302-2350 


HST/GHRS G270M 


1.1 


40 


08/24/92 


2350-2775 


IUE SWP-LO 


0.75 


70 


many c 


2775-2822 


HST/GHRS G270M 


1.1 


120 


08/24/92 


2900-5700 


HST/STIS G430L 


1.0 


100 


07/12/00 



"Scaling factor to normalize flux to the HST fluxes. 
b Approximate S/N in a 5 A bin. 

c Multiple IUE spectra were coadded by (Valenti, Fallon, Sz Johns-Krull 2003) and 
are available at http://archive.stsci.edu/prepds/iuepms/datalist.html. 
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Table 5: Properties of emission lines in the STIS spectrum" 



ID 


*obs 


V 


r WlilVl 


r 1 lux 


ID 


*obs 


V 


T7TTrTTT\ /T 


r lux 




X 

A 


km s _1 


km s 1 






X 

A 


km s 1 


km s 1 




("\ TTT 

L 111 


1 1 i /* i 
1175.164 




OOT 

327 


ctt i 

57.1 


T7\ , TT 

be 11 


1 o/^n i i o 

13b0.112 


1 o 

-12 


1 oo 

128 


n cr 

9.5 


/~1 TTT 
L III 


117b. 5bo 




i no 

103 


14.4 


^e II 


13bb.344 


-9 


i nT 

197 


1 o n 

12.0 


TT 
bl 11 


1 i nn £?/^o 

iiyu.obz 


a A 
64 


1 O A 

134 


o cr 

3.5 


T7V, , TT 

re 11 


i oro c cr O 

13b8.b53 


OO 

-32 


i no 

11)2 


o o 

3.2 


C " TT 
bl II 


1194.462 


o 

-8 


tct 

75 


O T 

3.7 


TT 1 TTO 

be II: 


1 OTn AT A 

1379.474 


o 

3 


01 n 

210 


o n 

8.0 


O • TTT fl \ d 
bl III (1) 


1206.149 


o r 

-85 


i nr 

105 


o < o 

24.8 


bi IV (1) 


i ono or'C\ 

1393.860 


o cr 

25 


o i n 

310 


o A n 

240 


O" TTT /0\tZ 

bi III (2) 


120b. 495 


1 


O/^O 

362 


O 1 /I 

214 


C " nr / o \ d 

bi IV (2) 


1 ono or* 1 

1393.861 


o cr 

25 


i r o 

158 


i m 

150 


Jri 1 


101^7 ran 

1217. 5UU 






75 


c : tt r / 1 ^ d 
bi IV (1) 


1 A no coc 

1402. b3b 


07 

-27 


1 1 O 

113 


70 

78 


IN V (1) 


i ooo i cr i 

12Jo.l51 


-160 


3b8 


oo n 

28.0 


c : tt r fo\d 
bi IV (2J 


1 A no 7fi A 

14U2.7U4 


1 O 

-12 


OTCT 

275 


oo cr 

225 


JN V (2) 


1 oo o ooo 

1238.828 


4 


i o cr 

135 


on o 

20.3 


T ^ t O 

F I: 


1/111 TOO 

1411.783 


1 o 

12 


111 

111 


r t 

5.7 


o 

f 


1 o r o i cr /i 

1252.154 




O 1 

81 


o n 

3.0 


TT 1 ttO 

l*e II: 


i/ii/i non 

1414.920 


1 n 

10 


1 TO 

172 


A cr 

4.5 


C; ttc 
bi 11 


12b5.000 






1 1 o 

11.3 


/-i TT 

L-r 11 


1 /lOO CO /I 

1428. b24 


T 

-7 


C A 

54 


o n 

2.0 


r \ T o 

o I: 


1280.95b 


O O 

28 


77 


i cr 

1.5 


T7\ , tt 

be 11 


1 AG A mo 

14b4.918 


1 1 
-11 


1 1 o 

118 


cr o 

5.3 


T7A TT 

J^e 11 


1 OOO d A 1 

128o.b41 


1 7 
-17 


OOO 

228 


cr o 

5.8 


CJ T 
b I 


1 /1 71 OAT 

1471.847 


r 

5 


i i cr 

115 


1 O T 

13.7 


C " TTT 
Si III 


1 on A A C 1 

1294.451 


on 

-20 


1 OT 

127 


i n 1 

10.1 


/"I TO 

LI; 


1 /I o i /i r 

1481. b45 


OO 

-22 


O 1 1 

211 


t r 

7.5 


C T 

b I 


1 on cr pi n 

1295. oil) 


-8 


i nn 

100 


oo o 

22.3 


C TO 

ai; 


1 /i o o /lor* 

1483.42b 


O 1 

21 


01 O 

213 


1 O O 

12.2 


C T 

b I 


129b. 05b 


-25 


oca 

2b0 


a a n 
4b. 9 


T7V, , TT 

be 11 


1 rrnt? C70 

150b. 572 


i n 

10 


1 07 

127 


O T 

8.7 


C; TTT 

bi III 


1298.992 


1 o 

13 


ono 

203 


oo o 

22.8 


c: ttc 
bl 11 


1 cot? nnn 

152b. 000 






1 /I o 

14.2 


/ \ tC 
U 1 


1 ono nnn 
1302.000 






91 


c; TTC 

bl 11 


1 coo nnn 
1533.000 






on 
2U 


(J I 


1 OH A T/"* o 

1304.762 


on 

-20 


A T 

47 


o o o 

33.3 


n tit / 1 \ d 

C IV (lj 


1 r a t nnn 

1547.900 


r a 

-54 


ono 

202 


1 TO 

173 


(J 1 


1306.000 






1 OT 

137 


U IV (2j 


1 cr /i o nno 

1548.092 


1 C 

-16 


A o cr 

435 


or O 

2b8 


c: ttc 

bl 11 


1 onn 1 

1309.261 


-1 


257 


en 

59 


ri Tir / n\d 
U IV (3) 


1548.287 


O 1 

21 


o c 

35 


Ol c 

31.5 


"D TT 

F 11 


1 onn C7n 

1309.870 


i 
1 


i n^; 

10b 


18.8 


U IV (1) 


i rrn onT 

1550.807 


9 


oon 

229 


1 no 

192 


C: ttc 

bi II 


1 01 n t at 
1310.747 


1 o 

12 


232 


1 c n 

lb.O 


ri Tir f rt\d 

Lj IV (2) 


i rrn nnn 

1550.839 


15 


b5 


A n 

40 


T~» TT r 

^e IP 


1 01 n on 

I3l2.b30 


O A 

34 


OOO 

238 


5.b 


T - ' TTT 

l*e II 


1 f f o no O 

1558.928 


OO 

-28 


lbl 


T O 

7.2 


C t 

b I 


1 oo o or cr 

1323.355 


O A 

-34 


1 A 1 

141 


A O 

4.2 


T -1 '. . TTC 

l<e II 


1 r o cr /i o 

15b3.543 


A r 

-45 


n 1 

91 


b.l 


C t 
b I 


1 one cr to 

132b. 572 


1 A 

-14 


no 

92 


i cr 

1.5 


TT 1 ^ TTC 

^e 11 


1 cr cr* a oc 

15bb.48b 


o cr 

35 


Ol 1 

211 


1 O O 

13.3 


P T 


lozy.oio 


1 1 
-ii 


1 80. 

ioy 


o.z 


TPo TT 
re 11 


10«jl.UU0 


1 7 
-1 ( 


1/1/1 


ion 
iz.y 


Fell? 


1330.985 


9 


14b 


4.8 


He II (l) d 


lb40.455 


-1 


25b 


35 


C II C 


1334.000 






b.3 


He II {2) d 


lb40.499 


7 


70 


28 


CII C 


1335.000 






112 


He II (3) d 


lb41.284 


150 


124 


12.7 


Cr II 


1347.13b 


24 


83 


2.9 


OIII] 


lb66.228 


lb 


237 


47 


CI I 


1351.b49 





135 


27.2 


AI IP 


lb70.000 






15 


01 


1355.bl5 


6 


157 


35.5 


Cr II? 


lb97.215 


-34 


184 


15.5 


CI? 


1357.240 


25 


93 


2.0 


Fe IP 


1702.000 






23.b 


01 


1358.438 


-14 


88 


5.5 













"Lines fit with single Gaussian, except where noted. fe 10 -15 erg cm~ 2 s _1 

c Wind absorption corrupts the emission line profile. rf Fit with multiple Gaussians 

e Possibly blended with Si I, Si II emission. labeled (1), (2), or (3) 
Possibly blended with S I 



-40- 



Table 6: Atomic and unidentified Emission Lines in the FUSE spectrum" 



ID 


Ao6s (A) 


v (km s 1 ) 


FWHM (km s" 1 ) 


Flux 6 


C III (l) c ' d 


976.154 


-263 


106 


3.7 


C III (2) c ' d 


976.679 


-102 


95 


10.5 


? 


997.313 




36 


0.8 


VI 


1031.251 


-190 


359 


30.4 


VP 


1036.861 


-215 


52 


2.8 


? 


1048.151 




67 


2.3 


? 


1074.385 




42 


1.3 


NI C 


1134.284 


-33 


128 


5.4 


NI C 


1134.974 





64 


3.3 


NI? 


1168.514 


-4 


79 


1.5 



"Lines fit with single Gaussian 
6 10~ 15 erg cm -2 s _1 



c Wind absorption corrupts emission profile. 
d Fit with multiple Gaussians. 
e Corrupted by H2 absorption. 
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Table 7: Pumped atomic lines in the STIS spectru m 

(A) 



ID 


Aofe (A) 




^pump 


Fe II 


1288.641 


HI 


1213. 


SI 


1295.610 


01 


1302. 


SI 


1296.056 


01 


1302. 


Sill 


1309.261 


01 


1304. 


PII 


1309.870 


01 


1301.874, 


P II 


1310.747 


01 


1305. 


Fe II 


1312.630 


HI 


1215. 


Fe II? 


1330.985 


I 


1304. 


Cr II 


1347.149 


HI 


1215. 


CI I 


1351.651 


CII 


1335. 


Fe II 


1360.117 


HI 


1215. 


Fe II 


1366.345 


HI 


1214. 


Fe II 


1368.649 


HI 


1214. 


Fe II? 


1379.474 


Si IV? 


1403. 


P I? 


1411.783 


HI 


1216. 


Fe II? 


1414.920 


cm 


1175. 


Cr II 


1428.624 


HI 


1217. 


Fe II 


1464.918 


HI 


1213. 


SI 


1471.840 


I 


1302. 


Fe II 


1506.544 


CII 


1335. 


Cr II? 


1697.215 


HI 


1304. 
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Table 8: Identification of H2 lines in the STIS spectrum 



ID 


^obs 


(A) 


Flux a 


ID 




(A) 


Flux a 


ID 


^obs 


(A) 


Flux a 


1-2 R(3) 


1202 


396 


3 


1 


1-4 R(6) 


1327 


506 


1.7 


1-7 R(3) 


1489 


497 


10.6 


1-5 R(3) 6 


1208 


787 


1 


8 


1-4 P(5) 


1329 


047 


-1.0 


4-8 P(8) 


1494 


917 


1.2 


0-2 P(3) 


1225 


469 


2 


1 


0-4 P(2) 


1338 


512 


13.0 


1-7 R(6) 


1500 


364 


8.6 


2-2 R(ll) 


1237 


511 


1 


2 


0-4 P(3) 


1342 


192 


12.9 


1-7 P(5) 


1504 


679 


18.5 


1-2 P(8) 


1237 


827 


4 


5 


4-6 R(3) 


1359 


031 


1.4 


4-9 R(3) 


1513 


489 


1.7 


4-4 R(3) 


1253 


625 


1 


5 


4-6 R(6) 


1371 


055 


0.7 


3-7 R(15) 


1513 


916 


1.5 


3-1 P(17) 


1254 


074 


1 


7 


4-4 R(17) 


1372 


014 


0.9 


blend 6 


1516 


122 


4.6 


1-3 R(3) 


1257 


780 


6 


4 


4-6 P(5) 


1372 


671 


1.3 


0-7 P(2) 


1521 


529 


1.8 


3-2 R(15) 


1265 


121 





7 


2-4 P(13) 


1379 


913 


1.1 


1-7 P(8) 


1524 


586 


8.8 


4-4 R(6) 


1266 


515 





6 


1-5 P(5) 


1387 


299 


3.4 


0-7 P(3) 


1525 


092 


1.5 


4-4 P(5) c 


1266 


823 


1 


7 


0-5 P(l) 


1396 


157 


0.9 


2-8 R(ll) 


1555 


810 


3.2 


3-1 P(18) 


1268 


806 





5 


0-5 P(2) 


1398 


883 


10.8 


1-8 R(6) 


1556 


796 


7.5 


3-4 P(l) 


1270 


527 





6 


2-5 R(ll) 


1399 


151 


1.8 


1-8 P(5) 


1562 


317 


6.9 


1-3 R(6) 


1270 


967 


3 


1 


0-5 P(4) 


1407 


154 


3.2 


4-10 P(5) 


1572 


548 


1.6 


2-2 P(13) 


1271 


128 


2 


4 


blend d 


1408 


807 


1.8 


1-8 P(8) 


1580 


573 


4.8 


1-3 P(5) 


1271 


869 


5 


1 


3-5 R(15) 


1418 


169 


1.1 


0-8 P(2) 


1582 


350 


3.1 


4-2 R(17) 


1273 


968 


1 


7 


1-6 R(3) 


1430 


945 


9.6 


2-8 P(13) 


1588 


709 


6.7 


0-3 R(0) 


1274 


495 


6 


4 


2-5 P(13) 


1434 


439 


2.2 


3-10 P(l) 


1591 


263 


1.3 


0-3 R(l) 


1274 


855 


4 


5 


1-6 R(6) 


1442 


794 


5.8 


3-9 R(15) 


1593 


191 


3.9 


0-3 R(2) 


1276 


222 


1 





1-6 P(5) 


1446 


051 


14.2 


2-9 R(ll) 


1602 


176 


4.9 


3-2 R(16) 


1277 


958 





6 


2-6 R(ll) 


1453 


024 


3.4 


4-11 R(3) 


1602 


543 


3.0 


0-3 P(2) 


1279 


407 


5 


9 


0-6 R(0) 


1454 


761 


3.7 


1-9 R(3) 


1603 


177 


2.5 


0-3 P(3) 


1283 


051 


6 


1 


0-6 R(l) 


1454 


894 


5.1 


4-11 R(6) 


1606 


202 


2.1 


0-3 P(4) 


1287 


632 





4 


0-6 R(2) 


1455 


847 


1.0 


1-9 R(6) 


1610 


884 


2.5 


4-4 P(8) 


1287 


886 





6 


0-6 P(l) 


1457 


348 


1.1 


4-11 P(5) 


1613 


645 


4.8 


2-3 R(ll) 


1290 


842 





8 


0-6 P(2) 


1460 


095 


7.6 


1-9 P(5) 


1617 


817 


2.3 


1-3 P(8) 


1293 


814 


4 


4 


0-6 P(3) 


1463 


712 


11.6 


3-9 P(17) 


1622 


084 


1.5 


1-4 R(3) 


1314 


570 


2 


6 


4-6 R(17) 


1465 


096 


2.1 


4-11 P(8) 


1622 


821 


3.8 


3-2 P(18) 


1320 


133 





7 


1-6 P(8) 


1467 


014 


6.2 


2-9 P(13) 


1632 


458 


3.1 


3-5 P(l) 


1324 


966 





8 


4-8 P(5) 


1476 


996 


1.0 


3-11 P(l) 


1636 


275 


3.9 



a 10 -i5 erg cm -2 s -i b Werner-band (C-X), pumped by O VI at 1031.862 A 

c Also includes emission from 4-1 P(19) d Blend of 1-5 P(8) and 3-4 P(17) 

e Blend of 0-1 R(0) and 0-1 R(l) 
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Table 9: Lyman and Werner-band H 2 lines in 


the FUSE spectrum 




ID a 


Kbs (A) 


v (km s 1 ) 


FWHM (km s" 1 ) 


Flux 6 


C-X 2-2 P(5) 


1058.255 


-22 


31 


1.3 


B-X 1-0 P(5) 


1109.273 


-8 


36 


0.4 


C-X 1-3 Q(3) 


1119.017 


-14 


71 


0.7 


C-X 0-2 Q(10)+ C-X 1-3 P(5) 


1127.096 




61 


1.5 


B-X 1-1 R(3) 


1148.632 


-16 


30 


1.4 


C-X 1-3 P(ll) 


1154.850 


-14 


36 


1.5 


B-X 1-1 P(5)+ B-X 1-1 R(6) c 


1161.774 




60 


3.4 


B-X 0-1 P(2) 


1166.176 


-18 


34 


0.9 


B-X 0-1 P(3) 


1169.687 


-15 


38 


0.8 


C-X 0-3 Q(10) + C-X 1-4 P(5) 


1171.916 




47 


3.4 


blend d 


1174.396 




122 


4.9 


C-X 2-6 R(3) 


1174.786 


-19 


25 


1.6 


C-X 1-3 Q(16) 


1177.208 


-14 


21 


0.6 



a B-X refers to Lyman-band, C-X refers to Werner-band 
6 10~ 15 erg cm -2 s _1 

C B-X 0-1 R(0) and B-X 0-1 R(l) may also contribute 

C C-X 1-5 Q(7), C-X 1-4 R(9), C-X 2-5 R(0), C-X 2-5 R(l), C-X 2-5 R(2) 
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Table 10: Properties of Wind Absorption Lines" 



ID 


Kbs (A) 


v]!, (km s J 


v max ( km 


s ) hi (cm ) 


if/ / & \ 
W\ (A) 


Si II 


1189.75 


-170 


-290 





-I OO 1 O f\ *~7 

1.20 ± 0.07 


AT T 

JN I 


1 i oo o o 

1199.83 









i r\ r* i a 1 1 P 

1.06 ± 0.11 


Si III 


1 OA/ 1 1 r 

1206.75 


i oo 

-180 


ooo 

-220 





O OO 1 O OO 

0.20 ± 0.02 


AT T 

JN I 


1242.74 


i i o 

-110 


i oo 

-190 


1 ooo A 

19224 


O CiC 1 O O A 

0.86 ± 0.04 


bl II 


1265 


-150 




OOT 

287 


1.15 ± 0.15-^ 


O I 


1 OA1 A O 

1301.43 


-170 


-290 





O O /I 1 O 1 

0.94 ± 0.01 


O " TT 

Si II 


-1 OAO f fl 

1303.57 


-180 


-260 





o 1 i o i r\ci 

0.61 ± 0.10 y 


I 


1304.30 


-130 


-220 


158 


0.72 ± 0.10 9 


(J I 


1 onr on 

1305. 3D 


-170 


-280 


227 


c\ oo i o m 

0.88 ± 0.01 


O " TT 

Si II 


1308.55 


-170 


-270 


OOT 

287 


o to i o or 

0.78 ± 0.05 


C II 


1334.5 


-140 


-230 





O O O 1 OOP" 

0.98 ± 0.25 


1 TT 

C II 


1 O O A O O 

1334.93 


-170 


-300 


63 


1 O O I O 1 

1.08 ± 0.01 


CI I 


1 o A or 

1346.35 


-200 


-320 





O O O 1 o ■ o 

0.38 ± 0.10 


O T 

S I 


1 A O I O 

1424.60 


-120 


-210 





O C\1 1 O O A 

0.27 ± 0.04 


AT T 

JN I 


1 /( oo oo 

1492.02 


-120 


O 1 o 

-210 


1 ooo J 

19224 


o 11 i o or 

0.77 ± 0.05 


AT T 

JN I 


1 A c\ o Tr 

1493.75 


i oo 

-190 


O 1 O 

-310 


1 ooo ^ 

19224 


o r c i o oo 

0.56 ± 0.03 


O " TT 

Si II 


1525.84 


-170 


-280 





o o i o o r 

0.96 ± 0.05 


O " TT 

Si II 


1532.55 


-170 


-280 


OOT 

287 


o o i i o o r 

0.97 ± 0.05 


Fe II 


1558.18 


-170 


-270 


1872 


o r* o i o oo 

0.60 ± 0.09 


TT 

re 11 


i no 
1562.93 


-16U 


-26U 


2431) 


U.5U ± U.U5 


Fell 


1565.95 


-70 


-140 


1872 


0.40 ±0.10 


Fe II 


1574.23 


-150 


-200 


1872 


0.23 ± 0.05 


Fe II 


1579 59 


-190 


-290 


2430 


0.53 ± 0.04 


Fe II 


1607.55 


-170 


-270 





0.66 ± 0.08 


Fe II 


1612.79 


-190 


-270 


1872 


0.55 ± 0.09 


Fe II 


1620.73 


-170 


-230 


385 


0.34 ± 0.04 


Fe II 


1624.64 


-160 


-260 


2430 


0.47 ± 0.06 


Fe II 


1628.15 


-180 


-250 


667 


0.28 ± 0.07 


Al II 


1669.85 


-170 


-270 





1.02 ±0.06 


Fe II 


1701.22 


-140 


-200 


1872 


0.53 ± 0.05 


"Measured assuming an emission line symmetric about the radial velocity of the star 
or as absorption against the continuum. 

^Central velocity of absorption. c Maximum velocity of absorption. 
^Energy of lower level. e Multiplet. Uncertain blue edge ff Blend 



